A little inflation at the cosmological QCD phase transition 
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We reexamine the recently proposed "little inflation" scenario that allows for a strong first order 
phase-transition of QCD at non-negligible baryon number in the early universe and its possible 
observable consequences. The scenario is based on the assumptions of a strong mechanism for 
baryogenesis and a quasistable QCD-medium state which triggers a short inflationary period of 
inflation diluting the baryon asymmetry to the value observed today. The cosmological implications 
are reexamined, namely effects on primordial density fluctuations up to dark matter mass scales 
of M max ~ IMq, change in the spectral slope up to M max ~ 10 6 Mq, production of seeds for the 
present galactic and extragalactic magnetic fields and a gravitational wave spectrum with a peak 
frequency around v pea k ~ 4 ■ 10 _8 Hz. We discuss the issue of nucleation in more detail and employ 
a chiral effective model of QCD to study the impact on small scale structure formation. 

PACS numbers: 



I. INTRODUCTION 



At about 10 microseconds after the big bang a phase 
transition from the quark-gluon plasma to a hadron gas 
is expected to have taken place at a temperature of about 
Tqcd ~ 150 — 200 MeV. In the last decade it has become 
more and more clear that this transition was most proba- 
bly only a rapid crossover as indicated by more and more 
refined lattice gauge theory calculations at zero baryon 
density [10]. In standard cosmology the baryon asym- 
metry is tiny t\b = tib/s ~ 10~ 9 , with tib being the 
net baryon density and s the entropy density, as deduced 
from later stages in the evolution of the universe. There- 
fore a first order QCD phase transition seemed very un- 
likely given the conditions. Still, the QCD phase dia- 
gram is for most parts terra incognita. The chiral and 
the deconfmcment transition do not necessarily coincide 
but there are some indications from effective models Q 
and lattice QCD calculations that there is at least a sig- 
nificant connection between the two. There has been 
recent progress in the attempt to include a finite baryon 
density on the lattice [1, [j| but effective models are still 
the method of choice to explore the uncharted regions of 
the QCD phase diagram Q. Findings indicate that at fi- 
nite baryon densities a first order phase transition can be 
expected as shown by chiral effective models of QCD Q 
caused by the melting of quark and/or gluon condensates 
or by color superconductivity Q. A sketch of a possible 
QCD phase diagram is depicted in figure Q] along with 
the commonly accepted path the universe took during 
and after the QCD-transition. The universe starts out 
in the upper left and move along the temperature axis 
from the chirally symmetric quark gluon plasma through 
a crossover transition to the chirally broken hadron gas 
phase. Once protons and anti-protons stop to annihi- 
late below 35 MeV the baryon chemical potential quickly 
shoots up from 1 eV to the nucleon mass (see ref. @ for 
more details). Effective models of QCD O El as well 
as lattice calculations Q at finite baryon chemical poten- 



tial give hints for the existence of a critical endpoint at 
M c = 0(1)T C . 
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FIG. 1: Sketch of a possible QCD-phase diagram with the 
commonly accepted standard evolution path of the universe 
as calculated e.g. in [i| depicted by the green path. 



At this point one might ask if there is a simple scenario 
with the cosmological QCD phase transition being first 
order without violating the constraint of a small baryon 
asymmetry in the later evolution of the universe. In a 
recent publication we have introduced the little inflation 
scenario [HI], that allows for such a first order QCD phase 
transition in the early universe without being in contra- 
diction to present cosmological observations. In figure [5] 
we sketch the evolution path of the universe in the little 
inflation scenario. Here the universe starts out at a large 
baryon chemical potential and therefore crosses the first 
order phase transition line but stays in the deconfincd 
chirally symmetric phase. The universe is trapped in the 
wrong QCD vacuum state and undergoes a short period 
of inflation until the delayed phase transition takes place. 
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The released latent heat then causes a large entropy re- 
lease that dilutes the baryon asymmetry to the presently 
observed value. Afterwards the universe evolves along 
the standard path just as in figure [T] 
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FIG. 2: Sketch of a possible QCD-phase diagram with the 
evolution path of the universe in the little inflation scenario. 



The concept of a little inflation (or tepid inflation) at 
the QCD phase transition has been also introduced ear- 
lier by Kampfcr ct al. [l3l - [l6j and for a inflationary pe- 
riod of similar duration as discussed here later by Borgh- 
ini et al. [I?} ■ In both cases an initially higher net baryon 
density is diluted to the presently observed small value 
of the baryon-to-photon ratio in the course of the QCD 
phase transition. The general idea of a short inflation to 
reduce a too high baryon asymmetry was mentioned even 
earlier by Linde in a publication on Affleck-Dine baryon- 
gensis [I8| but not explicitly in the context of the QCD 
phase transition. 

In the present manuscript we will reexamine the find- 
ings presented there in more detail and employ a chiral 
effective model of QCD to study some of the implications 
more thoroughly. The structure of the paper is as follows. 
In sec. II we address the first prerequisite of the scenario, 
namely a baryon asymmetry of order one before the QCD 
phase transition. In sec. Ill we discuss the second prereq- 
uisite, i.e. the topic of nucleation and especially the issues 
of supercooling and high surface tension. In sec. IV we 
introduce the dilaton quark meson model and apply its 
results in sec. V where we discuss the implications for lin- 
ear structure formation. Sec. VI deals with the changes 
to dark matter physics. Generation and modification of 
magnetic fields and gravitational waves are discussed in 
sec. VII and VIII, respectively. In the appendix we sum- 
marize some analytic limits for the structure formation 
calculations that are mostly not found in the literature 
to our knowledge and may help with the understanding 
of the numerical results in sec. V. 



II. BARYON ASYMMETRY 

As we have seen one of the main requirements of such a 
short inflationary period at the QCD phase transition is 
a non-vanishing baryochemical potential (J-b/T <~ 0(1). 
Big bang nucleosynthesis calculations predict the ob- 
served primordial abundances of elements correctly only 
if the baryon asymmetry was tiny at a temperature of 1 
MeV and below. The cosmic microwave background radi- 
ation as well as large scale structure observations predict 
very similar values and combining all these observations 
one finds a baryon asymmetry of 5.9 • I0~ 10 < rjs < 
6.4 • 10" 10 at 98% confidence 

Now we need to estimate how long such a little infla- 
tion has to be in order to start out with a sufficiently 
large ratio of fis/T. The net number of baryons in a 
comoving volume is conserved and can be estimated by 
Nb ~ o^nsiTf ~ Oj/iB/T 1 / where the index % refers to 
the initial values when the vacuum energy starts to dom- 
inate the energy budget of the universe and / to the final 
values after reheating. Therefore the initial ratio of the 
chemical potential to the temperature can be higher by 

VBi „ 3 HBf f T f\ 3 m 

with 9 = ctf/cii. If the phase transition at the end of 
inflation transpires on a timescale much shorter than the 
Hubble time then the universe reheats back to the initial 
temperature at the start of inflation in good approxima- 
tion Ti ~ Tf. Then we can conclude from equation ([!]) 
that for 8 ~ 10 3 « e 7 the baryon asymmetry before in- 
flation rjBi and jUj/Tj will be of order unity. The latter 
condition would, as we have seen, suffice to allow the 
QCD phase transition to be first order. 

The next question we need to address is if such a high 
initial baryon asymmetry is possible within one of the 
established baryogenesis mechanisms. Baryogenesis has 
been a long standing problem in cosmology ever since 
the pioneering publication by Sakharov [20[ and is still 
a very active field of research especially since it became 
clear that successful baryogenesis requires physics beyond 
the standard model, see |2fl425| for some extensive review 
articles of the field. 

In the well established Affleck-Dine mechanism of baryo- 
genesis [IH, [26| a large baryon asymmetry is much more 
natural than in other classes of models like for example 
baryogenesis via leptogenesis. In short the idea is that 
baryon- and lepton-number carrying scalar fields with 
very flat potentials can locally aquire very large expec- 
tation values. The Affleck-Dine mechanism can readily 
be incorporated into supersymmetric models [27j , where 
squark- and slepton-fields play the role of the baryonic 
scalar fields. Once supersymmetry is broken the flat di- 
rections are lifted and the scalar-condensates decay to 
standard model particles leaving a finite baryon and lep- 
ton asymmetry. In simple realisations the Afflcck-Dinc 
mechanism can easily produce a too high baryon asym- 
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metry for the standard cosmological scenario, thus ei- 
ther models with multiple fields or more sophisticated 
coupling terms have to be introduced to limit the ini- 
tial baryon number production or a subsequent reduc- 
tion is necessary. The latter could be achieved, as men- 
tioned earlier, by a large entropy release that dilutes the 
baryon to photon ratio to the right value observed today 
for example by an inflationary period (see e.g. ref. [H|). 
That being said Affleck-Dine baryogenesis can provide 
7]b ~ 0(1), where this is probably an upper limit fl8j . 
Still, this bound has to our knowledge not been explored 
any further after the estimates in the initial publications 
by Affleck, Dine and Linde for the obvious reason that 
an even higher baryon asymmetry was not desireable. 

Now we want to make the above rough guess for the 
highest possible [1b before such a little inflationary period 
a bit more quantitative taking rjs = 1 as an upper limit. 
Due to asymtotic freedom QCD should be well described 
by a free gas of quarks and gluons at sufficiently high 
energy densities so we can use this to gain some more 
quantitative estimates. To keep things simple we take all 
particles to be masslcss. The energy density, pressure, 
entropy density and number density of a relativistic gas 
read 
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Here g is the effective number of bosonic helicity states, 
i.e. fermionic helicity states are weighted with a factor of 
|. For ngwe can directly use equation Q with g = g q /3 
for the degrees of freedom. The entropy density has con- 
tributions from particles with sizable chemical potential 
and from those without, therefore we label the quark 
degrees of freedom an index q and those that have a non- 
negligible chemical potential with an index [i. This is 
necessary because both are not necessarily the same since 
leptons should most likely carry an asymmetry similar to 
the baryonic one. 



s = g—T 3 + g u — u, 2 T 
y 45 ^ 21/ 



(0) 



If we now combine both we arrive at an estimate for the 
baryon asymmetry 



2£g 



i ^ 



+ g^T gl^ 



(7) 



Interestingly this means that in the limit of /j > T as 
well as in the limit /i <C T the baryon asymmetry is just 




proportional to fi/T. The limits can be directly read of 
to be 



(8) 



Here we assume for simplicity that all particle species 
with a non-zero chemical potential have the same chemi- 
cal potential i.e. [i = /i g = /i„ = /i e et cetera. In the end 
we assume that the equilibrium condition for the quark 
and baryon chemical potential holds [Lb — 3/i g . Note 
that one cannot treat baryons as fundamental degrees 
of freedom satisfying equation (j4} with a charge of 1 /3 
within this simple estimate or there would be a contra- 
diction to the chemical equilibrium condition. 

In figure ([3]) the results from ([7]) are shown for two par- 
ticle compositions each for negligible lepton asymmetry 
and for equal baryon and lepton asymmetry. One can 
see that the influence from the particle composition is 
only a small effect, since the additional degrees of free- 
dom contribute in a similar magnitude to numerator and 
denominator of t\b- On the other hand t\b is significantly 
suppressed at the same chemical potential when adding 
an equal lepton asymmetry. This can be easily under- 
stood since the lepton asymmetry only increases the en- 
tropy but not the baryon number. The limiting values 
for pb/T assuming T)b = 1 for the four cases are shown 
in table [Q 
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for TjB = 1 


A 


47.75 


21 


21 


88.89 


B 


61.75 


31.5 


31.5 


88.74 


C 


47.75 


21 


29.75 


125.7 


D 


61.75 


31.5 


43.75 


123.1 



TABLE I: Degrees of freedom in the 4 considered cases A-D 
correspond to the curves in figure [3] 

Next we can translate the limits on the initial chemi- 
cal potential to temperature ratio to a constraint on the 
extent of the dilution. The baryon number in a comov- 
ing volume is conserved, i.e. ubi = 9 3 n,Bf, therefore the 
extend of the dilution can be directly inferred from the 
ratio of baryon asymmetries before and after inflation 
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Note that this definition does not necessarily coincide 
with the length of the period of exponential expansion 
defined by negative total pressure as we shall see later. 
To evaluate this expression we only need to calculate the 
baryon asymmetry rj B i because the two specific entropy 
densities Si and s/ are by definition equal. 
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FIG. 3: Here we plot the resulting baryon asymmetry from 
equation (0 as a function of hb/T for different cases. All 
curves include photons, three neutrino families, electrons and 
positrons, up and down quarks and eight gluons. Solid red 
(A) and dotted blue (B) curves assume a negligible lepton 
asymmetry while the latter also includes strange quarks and 
muons. The dashed black (C) and the dashed-dotted orange 
(D) lines include a lepton asymmetry and again the latter 
adds s-quarks and muons. 



If we now make use of the experimental value for rjBf 
we find the upper limit on the inflation length is given 
by 



1176 r, 1 ^ 



(10) 



independent of particle composition. In figure 2] we show 
the corresponding maximum dilution of baryon number 
by a delayed QCD phase transition in the little inflation 
scenario. This figure is part of the results of the dilaton 
quark meson model and the structure formation calcu- 
lation in sections IV and V, respectively, but it is quite 
model independent apart from the value of the chosen 
value of the vacuum energy. 

The period of exponential expansion could also be es- 
timated by comparision with the condition py + Pr = 0, 
i.e. the point at which the pressure turns negative. Even 
such a simple estimate turns rather lengthy and would 
also not be very accurate because for interesting infla- 
tion lengths the dark matter energy density becomes of 
similar order than the vacuum and radiation energy den- 
sities. We will show the numerical results for the period 
of exponential expansion in contrast to the length of di- 
lution as shown in figure [3] in section V. 



III. NUCLEATION 

The next critical requirement of the little inflation sce- 
nario is a large supercooling or in other words if a suffi- 
ciently delayed phase transition is possible. This issue is 
directly connected to the stability and height of the bar- 
rier between the chirally broken phase and the chirally 
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FIG. 4: The reduction of r/s to the presently observed value 
from an initial value of one. The scale parameter a is normal- 
ized to the scale parameter at which reheating occurs. 



restored phase in the effective potential for sufficiently 
low temperatures. In chiral models of QCD including 
gluonic degrees of freedom in the form of a dilaton field 
the barrier only vanishes in the T — > limit 28[ thus 
strong supercooling is in principle possible and wc will 
come back to this model later on. 

First let us consider the nucleation rate T of the low 
temperature phase inside the high temperature phase 



r = r Q e~ AF */ T 



(ii) 



where the functional form is that of a thermally acti- 
vated process as found by Langer in the 60s and 70s, 
e.g. [29(. To is in general a temperature dependent dy- 
namical prefactor and AF* is the free energy needed to 
produce a critical sized bubble of the new phase inside 
the old phase. What is meant by a critical sized bub- 
ble in this context? If the temperature is smaller than 
the critical temperature of the phase transition T < T c 
the system becomes metastable and statistical fluctua- 
tions produce bubbles of the low temperature phase with 
a radius R and a free energy of 



AF=^-( PH (T) 



PL (T))R 3 + 4irR 2 a s 



(12) 



Here ph(T) and Pl(T) is the pressure in the high and the 
low temperature phase, respectively, and erg is the surface 
tension. The first term describes the energy gained by 
transforming a spherical volume of radius R to the new 
phase while the second term gives the energy it costs 
to create the surface interface around the bubble. Since 
Pl{T) > ph{T) both terms have opposite sign and there 
is a critical radius i?* at which AF has a minimum 



2a s 



Pl{T) -Ph(T) 



(13) 



only bubbles larger than R* can grow, for smaller ones it 
is energetically more favourable to shrink and disappear. 
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One might just estimate Tq by T 4 for dimensional reasons 
but Csernai and Kapusta [301 ] found T in an effective 
field theory to be 



16 /cts\ 3 / 2 <jsVhR* 
~ 3tt VST J ^{Awy 



(14) 



which can easily be a few orders of magnitude smaller 
than the naive estimate. Here r\H and £h are the shear 
viscosity and the correlation length in the high T phase, 
respectively, and Aw is the difference in enthalpy density 
w = p + p between the two phases. 

The important ratio for the cosmological QCD phase 
transition is T/H, i.e. the rate of nucleation to the Hub- 
ble parameter. Once this ratio exceeds unity bubbles are 
produced abundantly and coalesce until the transition is 
complete. If T/H does not exceed one then bubbles of 
the low temperature phase will form and grow but the 
distance between bubbles increases so fast that the vol- 
ume fraction of the new phase stays small. 

We will in the following compare to work done by Cser- 
nai and Kapusta (3(| El for the QCD phase transition 
within the bag model to find if the nucleation rate can be 
sufficiently small compared to the Hubble parameter such 
that the phase transition will initially fail. In ref. [3l[ 
the authors found that the transition is completed very 
quickly with only marginal supercooling of about 1% be- 
low the critical temperature. In fact this result depends 
strongly on the value of the surface tension as which they 
took to be ~ 50MeV/fm 2 . This number originates from 
an older work of Kajantie ct. al [32j who calculated the 
surface tension at critical temperature and zero density, 
for which the transition is found to be a crossover by all 
recent lattice calculations. As one can see from equation 
([TTjl T depends exponentially on the value of the surface 
tension as well as on the free energy difference between 
both phases and especially the former quantity is in prin- 
ciple unknown at non-zero baryon density. 
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FIG. 5: Nucleation rate over the Hubble parameter for the 
lowest value of the surface tension for which the phase tran- 
sition would initially fail. 

Using the bag model with T c = 170MeV and the same 
parameters as in [3(| we looked for the lowest surface 



tension at which T/H does not exceed unity at least until 
its maximum at around ~ T c /2. This might already be 
overstreching the applicability of (TrTj) but it should still 
give a reasonable estimate of the surface tension needed 
for nucleation to fail. We find that the surface tension 
must indeed be very large and exceed 448 MeV/fm 2 ~ 
3.7 T c 3 using their high value of the bag constant of B = 
(235McV) 4 . If we however go to the lower end of values 
found in the literature, i.e. the original number B = 
(145MeV) found by the MIT group to fit hadron masses 
|33j , we find that a significantly lower as = 124 MeV/fm 2 
suffices. The resulting T/H in that case is shown in figure 
[5] The surface tension for the QCD phase transition 
at non-zero baryon densities can only be estimated by 
effective models since lattice gauge theor y ca lculations for 
this case are still in its infancies. In ref. [34| a reasonable 
range of as = 50 — 150 MeV/fm 2 is discussed but even 
smaller or larger values are not excluded in principle. In 
Ref. [HJ the surface tension was computed within the 
linear sigma model to be as low as 5-15 MeV/fm 2 . If one 
considers very high densities the surface tension for the 
transition from color superconducting phases to nuclear 
matter could reach values of 300 MeV/fm 2 (36|. In figure 
[6] the minimal surface tension needed for nucleation to 
fail is shown for the commonly discussed range of the 
bag constant. 
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FIG. 6: Minimum value of the surface tension as a function 
of the Bag constant at which T/H does not exceed unity. 

Our estimate only covers the initial failure to nucle- 
ate but it is clear that the phase transition has to occur 
after some limited supercooling (compared to ordinary 
inflation) of only about 7 e-foldings at most as we have 
seen. We stress that one should not take this estimate 
too far because both B and as have to be temperature 
dependent in general since both will in a field theoretical 
approach originate from the relative height and the shape 
of the barrier between the two minima in the effective po- 
tential. Finally T/H must exceed unity for inflation to 
end and the phase transition to proceed, for which the 
surface tension has to drop sufficcntly fast such that fluc- 
tuations can easily overcome the barrier. Another possib- 
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lity would be the complete vanishing of the barrier and a 
spinodal decomposition as studied for example in (l5j for 
a bag like model. Other authors have also discussed the 
strong sensitivity of nucleation rates for example in the 
context of neutron stars and core-collapse supernovae. 
There it was found that nucleation timescales can basi- 
cally not be constrained and range from /is up to the age 
of the universe 



37] 



Also for heavy ion collisions strong supercooling is dis- 
cussed for the "quench" -scenario, see e.g. (38|. There the 
chiral phase transition is delayed as the field is trapped 
in a mctastable minimum and is only released to the true 
minimum in the T=0 limit. 

The equation of state has to fulfill the usual condition 
p + Sp < to enter an inflationary phase. In the bag 
model this would be the case below a temperature Ti n f = 

(30B/ (gn 2 )Y^ 4 . In the linear-cr-model or the NJL-model 
this occurs when the thermal contributions to the energy 
density become smaller than the vacuum contributions 
like the quark condensate (m q qq) w f 2 fn 2 and the gluon 
condensate /3 QCD /(2g) (G^Gf?) « 4B. 

We can conclude that QCD at non-zero baryon densi- 
ties is only poorly constrained and a delayed chiral phase 
transition is very well possible and has already been dis- 
cussed for several other scenarios apart from the early 
universe. 



IV. DILATON QUARK MESON MODEL 

To describe the dynamics of the phase transition and 
especially the impact on density perturbations it is es- 
sential to have a reasonable thermodynamic description 
of the chirally restored quark phase. For this we use 
the quark meson model with a dilaton field, which in- 
corporates chiral symmetry breaking as well as the trace 
anomaly of QCD. This model has been discussed by nu- 
merous authors [39l - l42j to describe nuclear matter. It has 
the interesting property that for a wide range of param- 
eters the high temperature phase does only disappear in 
the T — >• limit [28|, which is necessary to get an equa- 
tion of state of the "wrong vacuum" that can be used to 
model inflation. We will apply a simplified version of the 
lagrangian used in [4l| and stick closely to their nota- 
tion, while we do not include the w-meson for simplicity 
and use quarks instead of nucleons as degrees of free- 
dom. The La gran gian includes the linear cr-model first 
introduced by |43| which incorporates the scalar isovector 
7r-ficld with the cr-field as its chiral partner. Furthermore 
we include the isoscalar dilaton field x that incorporates 
the scale anomaly and thus a non-trivial vacuum of QCD. 
The Lagrangian reads 



2 ( d vxf 



where the potential is given by 
[/( 7 r,a, X ) = ^(a 2 + 7 r 2 ) 2 
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(16) 



that 



Note that the choice of the term —f^m^a 
breaks chiral symmetry explicitly is not unambigous. 
The given choice X 2 /Xo corresponds to a fermion mass 
term, while x/xo would lead to a bosonic mass term. 
The differences caused by this choice are small since the 
symmetry breaking term logarithmic in x is usually much 
larger, but for a similar model including vector mesons 
the quadratic choice is favored when comparing to nu- 
clear matter (40j . 

After integrating out the quark degrees of freedom [44[ 
one arrives at an effective mesonic lagrangian 



L(n,(T,x) = 2 ( 9 M 7r ) 2 + \ OV) 2 + \ (^X) 2 

U{k,<t,x) ~ ^qq(T,p,,m q ) (17) 



with the quark-antiquark potential reading 
n qq (T,p,a) = ^ ^dpp 2 [\n(l + e-^-^) 



o 



+ In (l + c-^+f) 



(18) 



Here the single particle energy is as usual given by E q 



p 2 + m 2 . The effective quark mass is on the mean field 

level determined by m 2 = g 2 a 2 . The glueball and meson 
fields have the following mean values in the vacuum (x) = 
Xo, (f) = do = fir and (ir) = 0. The full thermodynamic 
potential is then given by 

n(T,fi,-K,a,x) = (C/(tt,ct,x) - U vac )V + n qq (T,p,cr) 

(19) 

Where U vac is subtracted to ensure the correct normal- 
ization fl(0, 0, 0, f n , xo) = 0, yielding 



X ,4 ko 2 r2 2 i 7 



(20) 



i> W - 9(0 + il5T- tt)} V> - U(tt, a, x) (15) 



VL/V exhibits two minima, one at a ~ and x ~ 0.8xo 
and a second one at a ~ f„ and x ~ Xo- The former 
corresponds to the chirally restored phase with a low ef- 
fective mass while the second one is the chirally broken 
phase with a large effective mass. Full restoration of 
scale symmetry, i.e. the first minimum being located at 
a ~ and x ~ 0, is only realized at high temperatures 
and low densities at least for the parameters considered 
in our calculation. If the scalar coupling is sufficiently 
large both minima are present in the low temperature 
limit, although the chirally brokenphase is energetically 
favoured. As the authors of Ref. [42| have found the chi- 
rally restored phase will undergo a crossover to restored 
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scale symmetry at much higher temperatures if the den- 
sity is non-zero, i.e. the maximum moves towards x ~ 0- 
Note that we do not fix the effective quark mass in the 
vacuum via the Goldberger-Treiman relation because the 
model is set up to describe quarks in the high temper- 
ature chirally restored phase. We will later on fix the 
model parameters using the more relevant vacuum en- 
ergy constraints and the masses of the sigma meson and 
the dilaton. The constants kg and k\ are determined by 
the conditions 



an/v 



da 



dx 



= 



(21) 



The equations of motion are found by minimizing the 
thermodynamic potential with respect to a and %. 



dSl/V 
da 



on/v 

dx 



= = Act 1 - kn 



gps 



a - Uml 



i X 2 o r 2 

-k —a - 2/ 7r m 7r <7 2 



X 



Xo 
4fci 

Xo 



X_ 



1 + ln 



(22) 



(23) 



These equations can reduced to a one dimensional prob- 
lem by solving for x explicitly 



X = Xo 



Xa 3 
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k Q a + m% f w J 
here ps, the scalar density, is defined by 
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(24) 
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The pressure is as usual just given by 

the net quark density is calculated via 
dfl/V 



e P(E q +») + i 

(25) 
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The energy density is then given by 
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The entropy density can as usual be deduced from the 
Euler-equation 



e = Ts — pV + iin q 



P - pn q 



(29) 



Calculating the speed of sound is a bit more involved but 
straightforward. By definition the speed of sound is the 
isentropic derivative of the pressure with respect to the 
energy density 
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Isentropic means nothing else but 
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Using the total differentials of pressure and energy den- 
sity and equation pip we arrive at the speed of sound in 
the form 
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(32) 



The remaining parameters xo and A are fixed via the 
QCD vacuum energy and the mass of the sig ma meson. 
Investigating the QCD trace anomaly Ref. [45[ found that 
at tree level the trace anomaly of QCD is given by 



% = liQCD/{2g) (G^Gr) 



(33) 



where O^" is the energy momentum tensor, /3qcd is the 
beta-function of QCD, g is the strong coupling constant 
and G^„ is the gluon field strength tensor. The trace 
anomaly of QCD then relates the vacuum energy to the 
parameter xo [ill 
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(34) 



We neglect the contribution from the first term repre- 
senting the quark condensate for simplicity as done in 
[4l| because its contribution is much smaller than the one 
from the gluon condensate given by the second term for 
the parameters we will use later on. QCD sum rules sug- 



gest \e v 



(240 MeV) 4 (see ref. |46j) while bag model 



estimates range from (235 MeV) 4 down to (145 MeV) 4 
in the original paper of the MIT group [33| . This results 
in a possible range for the parameter xo of 205 MeV 
< Xo < 339 MeV. We choose a cr-mass of 642 MeV, a 
dilaton mass of 1.5 GeV and a vacuum energy of (236 
MeV) 4 to achieve a critical temperature of 170 MeV for 
the phase transition at zero net density. The scalar cou- 
pling is chosen to be g = 7.5 which is the limiting value 
above which the chirally restored phase is present even 
in the T —y limit approximately. 

In figure [7] we show the resulting speed of sound and 
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FIG. 7: Square of the isentropic speed of sound and equation 
of state for the most extreme case (Mb/T — 125.7. 
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FIG. 8: Pressure p as a function of the a— field in units of 
/,r at the phase transition temperature T = 10.1 MeV and 
jis/T = 125.7. 



equation of state in the maximum case p-Bi/Ti = 125.7 
As one could expect the speed of sound stays very close 
to the relativistic gas value of c 2 = 1/3 because the effec- 
tive quark mass stays low in the chirally restored phase. 
The equation of state nicely interpolates from a relativis- 
tic gas (w = 1/3) to that of vacuum energy (w = — 1). 
The small kink in the speed of sound is caused by the 
merging of a third always metastable intermediate phase 
with the chirally restored phase which causes a sudden 
but small change in the effective mass. The existence of 
this third maximum in the pressure within this model 
has been discussed before for example by Mishustinet 
al. [391 ]. it can also be seen in figure [8] at a ~ 0.5/ w . 
There we show the pressure as a function of the cr-field 
(or equivalently the effective mass) at the phase transi- 
tion temperature T = 10.1 MeV and pb/T = 125.7. The 
first minimum at a ~ is the chirally restored phase, 
the chirally broken phase is located at a ~ f n . The 
intermediate phase only appears close to the phase tran- 
sition and never becomes the favoured one. Note that at 
these temperatures and densities one may expect color- 
superconducting quark matter in one of many possible 
phases [|[ which exceeds the scope of the current inves- 
tigation but may be an interesting starting point for an 
alternative field theoretical description of the scenario. 

In the next section we will use the e,p, c 2 and w of the 
chially restored phase for our structure formation calcu- 
lations. 



V. STRUCTURE FORMATION 

Next we will investigate the effect of a little inflationary 
period on primordial density perturbations. In particu- 
lar dark matter perturbations are affected in several ways 
and on much larger scales than usual for the cosmological 
QCD phase transition. First of all the Hubble radius is 
roughly given by Rh ~ g~ x l 2 vnp{F~ 2 ~ 10 km which en- 
closes a total energy corresponding to about 1M Q . Since 



this epoch is long before matter radiation equality, i.e. 
Pdm ~ iflQCD / o-eq)pr ~ 10~ 8 p,R, the mass of dark 
matter in the same volume is smaller by the same factor 
resulting in a dark matter mass scale of approximatly 
10 -8 M Q . About ten years ago Schmidt, Schwarz and 
Widerin investigated the effect of the QCD phase transi- 
tion on dark matter perturbations [13, 13 . They found 
that peaks and dips in the spectrum of dark matter per- 
turbations may form for a first order phase transition but 
even for a crossover one could expect a boost for small 
scale perturbations. These effects were due to the re- 
duction of the speed of sound c s and equation of state 
w = p/p of the radiation fluid during the phase tran- 
sition. As the above estimate implies they only found 
these effects at very small mass scales below the Hubble 
scale. We shall examine the little inflation scenario with 
the same approach to density fluctuations, i.e. we work 
in the so called uniform-expansion gauge (UEG) that is 
free of spurious gauge modes and well behaved in the 
superhorizon limit [48l . |49| which is especially important 
for the little inflation scenario since modes can enter the 
horizon and exit the horizon several times. 

For the background evolution we assume a flat 
Fricdmann-Robcrtson- Walker metric that implies the 
well known Friedmann equations 



H = -47rG(p + p) (36) 

The overdot as usual denotes a derivative with respect 
to time. The common approach now is to decompose the 
perturbed metric into scalar, vector and tensor perturba- 
tions where only scalar perturbations lead to the growth 
of structure. The corresponding perturbed line element 
reads 

ds 2 = a 2 {if) [(1 + 2a) drj 2 - 2B\ i dx i dr] 

- ([1 + 2<p] 7 . y - + 2E\ ij ) dx l dx j ] (37) 
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Here a, B, ip, E are spacctime dependent scalar functions 
representing the four scalar degrees of freedom and 7,;,- is 
the spacial part of the background metric. Following [50| 
we introduce two combinations of the metric variables 
that are independent of spacial gauge transformations 
just like a and ip. These are 



X = -a(B - Ha 2 E') 
K — 3H (a + cup ) + - 



(38) 
(39) 



where primes denote derivatives with respect to the scale 
parameter a and k is a comoving wavenumber. One can 
show that x an d k describe the perturbations in the shear 
<Jij and the expansion O, respectively. The latter is given 
by = 3H — k, therefore choosing n = corresponds 
to having an unperturbed Hubble flow. This explains 
the naming uniform expansion gauge or Hubble constant 
gauge. Choosing x = leads to the more popular lon- 
gitudinal or conformal Newtonian gauge as used in the 
well known review of Mukhanov, Feldman and Brandcn- 
berger (Hlj . For ideal fluids the evolution equations in 
UEG read S M 



e = -3i?(e + 7f) - Ai;-3H(p + p)a (40) 
ip = -3Hip-TT- (p + p)a (41) 

which can be deduced from energy-momentum conserva- 
tion and the three divergence of the Euler equation. Here 
e = 5p and ff = Sp denote the perturbation of the energy 
density and pressure, respectively. Furthermore ip is the 
potential of the momentum density S, i.e. Vip = S. The 
latter is related to the fluid velocity v via ipk/a = (p+p)v. 
Equations (|4"D"jl and (j4"Tj) apply for each decoupled ideal 
fluid. All fluids are gravitationally linked via the pertur- 
bation of the lapse a and Einsteins i?Q-equation 



are connected via the last equation for the perturbation 
of the lapse a. The mean density contrast, equation of 
state and speed of sound arc calculated by 



5 = 



E t 5 iPi 

E* pi 



w 



Ei Pi 



E; hPi 



(46) 



Basically all viable dark matter candidates are already 
chemically decoupled from the radiation fluid at the QCD 
phase transition, thus their numbers are not repopulatcd 
by reheating after inflation. Therefore the dark matter 
number density is diluted by the same factor 3 as the 
net baryon number. As stated before the dark matter 
mass enclosed inside the Hubble horizon is of the order 
of l(r 8 Af at Tqcd ~ 170 MeV. Thus any influence 
on perturbations inside dark matter would not have any 
consequences on larger scales. An inflationary period at 
the QCD-phase transition can change this in two ways. 
First of all the amount of dark matter enclosed inside the 
horizon must be larger by a factor 9 3 initially to match 
the present day dark matter density despite the dilution. 
For a short inflationary period, as discussed here, one en- 
counters an additional effect on perturbations that have 
physical wavenumbers k p h < H at the beginning of in- 
flation. This is caused by an additional scale apart from 
H, namely if 1 / 2 , via equation (|42|) emerges. One may 
realize this by combining equations (|35[) and (|36[) to find 
that 



H 
IP 



It 
3 



P 



P 



w) 



(47) 



So as long as w is not too close to -1 both scales coincide, 
but during an inflationary phase this is no longer true. 
Let us do an estimate for a general mix of radiation, dark 
matter and vacuum energy. In this case 



(A + 3H)a = 4ttG(p + 3p) 



(42) 



Here one already realizes that H 2 and H appear as two 
distinct scales in the set of perturbation equations. These 
will be similar for most cases, but during an inflationary 
period they are not as we will see later. Introducing 
dimensionlcss variables 6 — Sp/p, tp = ktp/(ap) and the 
equation of state w = p/p the UEG set of equations takes 
the form 



6'i - 
# = 



3(c 



a 



^^ + —^-3(1 + ^- (43) 

a Ha a 

-i>i - cl — S, - (1 + w^ — a (44) 

a Ha Ha 

t( 1 + 3 ""> ,45, 



(1 



Here the index i refers to an individual fluid, each of 
which has a set of equations (|43| and (|44[) . H = Ha is the 
so called conformal Hubble parameter. Primes denote 
derivatives with respect to the scale parameter. All fluids 



H = -AttG 



~Pm — 
6 V a 



'a l 

Pui\ — 
a 



where the subscripts refer to matter and radiation with 
q = 3 to 4, respectively. The index i refers to the onset of 
inflation. Comparing this to the first Friedmann equation 
one finds that 
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As a consequence the two scales differ by 
1/2 



H 



H 2 



a 



-Pv 
(49) 



(50) 



This would not play any role for a long inflationary pe- 
riod, i.e. with more than 50 e-foldings. In this case H^ 1 ^ 2 
is beyond the size of the observable universe, approx- 
imately at the order of the infrared cutoff of the pro- 
duced primordial spectrum. Summarizing, there should 
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be two distinct scales in the spectrum dividing it into 
three regimes 
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One could even expect another spectral region of modes 
that are always sub-Hubblc till the end of inflation. 
These would be located below 10 -8 Mq and are thus of 
little relevance for structure formation. 
Translating our previous estimates to the highest affected 
mass scale involved we find 



M„ 



10- 8 M Q 9 3q/2 ~ (10 5 - 1O 9 )M (51) 



at most for d tn ' ~ 640. In this case the second mass scale 
between the first and the second spectral region could be 
expected at 



M med ~ 1O" 8 M 9 3 ~ (1 - 1O)M 



(52) 



For the numerical treatment we assume a scale invariant 
primordial Harrison-Zeldovich spectrum to be present be- 
fore the phase transition. Each wavenumbcr k is fol- 
lowed separately from a point where is was sufficiently 
super-hubble to apply the initial conditions for a radi- 
ation dominated universe given by the growing super- 
horizon modes [M K3(l 



A A ( k 



1 A ( k 
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9 / 



(53) 
(54) 



For the description of the radiation background we take 
the input from the dilaton-quark-meson model as de- 
scribed in the previous section. We add a massless ideal 
gas of photons, gluons, e ± and three neutrino families as 
in case C described in section HU For the dark matter we 
assume a decoupled pressurclcss non-rclativistic gas with 
a vanishing speed of sound. 

In figure [9] we show the resulting spectrum of primordial 
fluctuations after a little inflation in comparison to the 
spectrum as expected without little inflation. The spec- 
trum is given in terms of the transfer functions defined 
in the following way 



T R (k) 
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S 2 n + i' 2 R /cl R 

\ !/ 2 
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(55) 
(56) 



The "in" quantities are evaluated at (final) horizon en- 
try in the limit of small wavenumbers [48j . i.e. in the 
unaffected part of the spectrum for the little inflation 
calculation. Note that a transfer function of unity does 
not correspond to non-linear perturbations. For a scale 
invariant spectrum it simply corresponds to a amplitude 
of 10~ 5 — 10 -4 since the perturbations in radiation are 
frozen until the decoupling of the cosmic microwave back- 
ground radiation (ignoring low scale damping effects). 
The fluctuations are evolved for 12 orders of magnitude 
in a to ensure that the whole spectrum is super-horizon at 
the start of the calculation and is completely sub-horizon 
at the end. The final temperature for the parameters 
used is T end ~ 150cV a T reheat /10 6 . 



All scales below M„ 



100 



10 Mq show a suppression, 
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FIG. 9: Spectrum of primordial fluctuations dependent on the 
dark matter mass scale in units of solar masses normalized to 
the amplitude at horizon entry. As input we used the equation 
of state of the chirally restored phase of a dilaton quark meson 
model with equal lepton and baryon chemical potentials (case 
C). For comparison the spectra without a little inflationary 
phase is shown with dashed lines. 

those below Mh ~ Mq show additional features de- 
pending on their phase during horizon exit. Above this 
scale the spectrum of density perturbations is given by 
the primordial spectrum of density perturbations, e.g. a 
nearly scale invariant spectrum. In the appendix we also 
give the approximate analytic solutions for the differ- 
ent regimes. Summarizing they show that in the un- 
affected region the modes grow similarly to the radiation 
dominated super-horizon solutions while the intermedi- 
ate modes are frozen. This consequently explains the 
relative suppression of the intermediate modes compared 
to the large scale limit. 

The numerical result for the maximum mass scale is quite 
close to the lower bound in aboves estimate (|5T|) because 
dark matter has to be more abundant than radiation dur- 
ing almost the complete duration of the little inflation, 
which can be seen in figure [TU] Still this mass scale is of 
cosmological interest as it is comparable to that of glob- 
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ular clusters (GC) which were the first objects to form 
during primordial galaxy formation (for a comprehensi- 
ble overview of the topic the reader may have a look 
at the review by Harris |52]). Globular clusters are very 
compact star clusters of several hundred thousand to sev- 
eral million stars, with a radius of only ~ 10 pc. They 
are very metal poor objects and age estimates from stel- 
lar evolution models strongly suggest that they should 
already have been created during the formation of their 
host galaxy. Their mass function has a well defined peak 
at M gc ~ 2 • 10 5 Mq in contrast to younger star clusters 
whose mass function shows a steep power law distribution 
between 10 4 M Q < M yc < 1O 7 M with an index of s» -2 
[53| . There have been attempts to explain the preferred 
mass scale for GC by a higher Jeans mass at low metal- 
licity that preferred more massive clusters at early times. 
Other explanations include disruptive processes (for low 
mass clusters) and mass loss due to stellar evolution (for 
high mass cluster) that might produce a preferred mass 
scale as seen in n-body simulations [54[ starting from a 
steep mass function like the one of present young clus- 
ters. The latter point was dismissed by Vesperini et 
al. [55j whose simulations have shown that only very fine 
tuned initial conditions will result in GC properties that 
fit the observations if one assumes GC formed just like 
clusters do today. On the other hand they showed that 
an initial mass function that was almost flat below the 
present peak mass scale succeeds in reproducing the ob- 
servational data. Interestingly a power law cutoff in the 
dark matter fluctuations below a mass scale 1O 5_6 M as 
found in the little inflation scenario could thus help to 
explain this standing problem in galactic astrophysics. 
The suppression of power spectrum below these scales 
could also be interesting to study because of its impact 
on the cuspy core density distribution of dark matter in 
small galaxies and the large number of halo structures 
seen in standard structure formation. 
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FIG. 10: Evolution of the background energy densities of radi- 
ation plus vacuum contributions, dark matter and of radiation 
alone. As one can see the inflation length is only Q m * ~ 640 
while the dilution factor is 9 — 1176 



VI. DARK MATTER 

Apart from the impact on small scale structure formation 
for dark matter the little inflation scenario has further 
direct consequences on properties of dark matter candi- 
dates. For cold dark matter the dilution of the energy and 
number densities leads to the possibility of a matter dom- 
inated phase before the inflationary phase since the dark 
matter energy density after reheating is basically fixed by 
the present day value. This can also be seen in figure [T0l 
that displays the evolution of the different contributions 
to the energy density, wherein dark matter just starts to 
dominate right before the onset of inflation. To account 
for the different ratio of radiation and baryon densities 
before little inflation the dark matter density has to be 
larger by the same factor 8 3 as the baryon density. For 
9 > 10 3 the dark matter contribution actually becomes 
larger than the vacuum contribution which would in any 
case limit the maximum length of the exponential expan- 
sion to 
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(57) 



where the bag constant B represents the vacuum 
contributions of QCD. Interestingly this limit and the 
previously discussed limit from the Affleck-Dine baryo- 
genesis coincide by chance, while the latter actually 
limits the entropy release the former only limits the 
length of exponential expansion. As a side remark, to 
produce a complete spectrum of primordial fluctuations 
one would require 9 > 10 10 , far beyond both limits so 
little inflation cannot replace standard inflation. 
Still figure [TU] shows that the period of exponential 
expansion 9 ln f sw 640 is even shorter than this estimate 
because the energy density of radiation increases so 
strongly with the baryon asymmetry at a fixed vacuum 
energy. This difference in 8 and 8 m ' is caused by the 
different dependencies of the entropy and the energy 
density on /i or rather uq. 

What does a larger dark matter density before the QCD 
scale mean for the properties of cold dark matter? For 
non-relativistic decoupling of dark matter the weak inter- 
action cross section will no longer give the right amount 
of dark matter today. This is due to the fact that the 
dark matter annihilation cross section has to be much 
smaller, i.e. 



annih 



T weak 



9 3 



because Qdm oc 



1 



dm 



where we ignore logarithmic dependencies on the dark 
matter mass. This allows more dark matter particles to 
survive annihilation before freeze-out and thus increases 
the CDM number density before little inflation. This 
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gives the interesting prospect that the little inflation can 
be probed by ongoing and future collider experiments like 
the LHC since the discovery of a standard weakly inter- 
acting massive particle as the neutralino would exclude 
the scenario. 

Another case would be thermally decoupled ultra- 
rclativistic particles where the dilution of dark matter 
number densities can be incorporated in the ordinary 
temperature relation to the radiation background. Here 
little inflation leads to an effective shift in the tempera- 
ture relation 



T = T, 



1/3 



Dec) 



(58) 



This in turn modifies the relation of warm dark matter 
relic mass and decoupling degrees of freedom to match 
the present day density found for example in (56j 
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(59) 

This shifts the suitable mass of a thermal relic particle to 
a much higher value without the need for a large number 
of additional effective degrees of freedom at decoupling 
beyond those of the standard model. 
There can also be effects for baryonic dark matter as dis- 
cussed by Jcdamzik [57j | . During a first order phase tran- 
sition the speed of sound vanishes and thus sufficiently 
nonlinear density fluctuations can collapse during that 
time. For an exponentially small fraction of Hubble vol- 
umes that are overdense enough primordial black holes 
(PBH) may form. The mass spectrum of these PBH will 
be strongly peaked around 1M© which corresponds to 
the total (not just the dark matter) energy density in- 
side the Hubble volume at the phase transition. The 
produced abundance of PBH depends on the spectral in- 
dex and amplitude of the density fluctuation spectrum, 
which we have seen is different and in general more com- 
plicated in the little inflation scenario. Nevertheless it 
seems quite clear that the suppression of small scale den- 
sity fluctuations will also strongly reduce the production 
of such primordial back holes during the phase transition 
at the end of a little inflation. 

During the nucleation process lumps of quark matter 
or small quark stars could be produced but only with 
A/ ~ 1O _9 M0 as we argue that nucleation starts after 
the little inflationary epoch. 



VII. MAGNETIC FIELDS 

A standing problem in astrophysics is the origin of 
large scale magnetic fields that have strengths of up to 
B^ bs = O.lfiG on extragalactic and up to 10 fiG on galac- 
tic scales. To understand the existence of such magnetic 
fields with correlation lengths of typically 0.1 Mpc it is 



necessary to have an initial seed field generated before or 
during galaxy formation. The required strength of such 
seed fields varies strongly with the assumed amplification 
mechanism and may vary over many orders of magni- 
tude 10~ 30 G < Bf ed < 10- 10 G, see |H| and references 
therein for an overview of the topic. The seed fields may 
be generated during ordinary inflation or at a first or- 
der phase transition. The latter has been discussed for 
the QCD phase transition by numerous authors [59l - [6l| 
at a time when the phase transition at small baryon- 
asymmetry was still believed to be first order. The es- 
tablished mechanism for magnetic field production was 
the collision of hadronic bubbles during the phase transi- 
tion [6(|. Different masses of quarks and nucleons would 
lead to a diffusion of baryon number via the bubble walls 
and consequentially a baryon contrast close to the phase 
boundary would develope [59l |60| . This baryon contrast 
can be estimated by the ratio of the net baryon numbers 
in the two phases to be 



R = — 5- 



(60) 



Because muons and strange quarks are already slightly 
suppressed at the critical temperature T c the baryon con- 
trast would also cause a charge dipole layer at the phase 
boundary to develope. The resulting net positive charge 
density is 



p c = e (2/3n u - l/3nd - l/3ra s ) = Pens 



(61) 



where the indices u, d and s refer to the different quark 
flavours and the factor of proportionlity j3 depends on the 
temperature, chemical potential and the masses of the 
particles. For a small t\b and reasonable strange quark 
and muon masses f3 ~ 10 -2 — 10~ 3 . After a strong su- 
percooling muons and strange quarks will be suppressed 
resulting in /3 w 0.2 for the little inflation case. Cheng et 
al. estimated the magnetic field generated by the collision 
the hadron gas bubbles to be 
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(62) 



due to turbulent charged flow. Here the flow paral- 
lel to the bubble walls was assumed to have veloci- 
ties v ~ r n HqcD giving the main contribution to the 
field generation. The thickness of the the baryon ex- 
cess layer r d was estimated according to results of [62| to 



be rd 



'diff 



jr n with rdiff being the baryon diffusion 



length and r n the mean separation of nucleation sites. 
R/3 should be at least 0.3 with the above estimates up to 
values of ~ 10 — 100 if baryon number can be effectively 
piled up by the expanding bubble walls. Thus we arrive 
at magnetic fields of strength Bqcd = 10 8 — 10 10 G for 
low baryon asymmetry, i.e. for the standart scenario as- 
suming a first order phase transition. If the baryon con- 
trast exceeds R ~ 10 then any initial field may be readily 
amplified by magneto-hydrodynamic (MHD) turbulence 
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to the equipartition value (see [6l[ and refs. therein) 

B eq = yW 4 ^ (63) 

where Vf is the fluid velocity. Now we shall modify these 
estimates for the little inflation scenario. First of all the 
initial value of the baryon number contrast R between the 
two phases can be much higher because quarks are much 
more favourable carriers of baryon number than nucleons 
at such low temperatures of T ~ 17OMcV/0 ~ 0.2MeV 
at the end of inflation. The diffusion length will also 
be larger because both baryon and antibaryon densities 
tib , ng will be additionally diluted by a factor 9 3 result- 
ing in 

rdiff oc \ j\jn B + n B ~ 4^m 6 3/2 ~ 10 cm (64) 

for a random walk approximation. Thus devclopement of 
MHD turbulence should be expected resulting in equipar- 
tition of the magnetic field with a strength of B eq w 
10 12 G. The fluid velocities were taken to be Vf ~ 1 be- 
cause the released latent heat is much larger than the 
thermal energy. 

Next one may ask if such a strong magnetic field does 
not violate bounds for the total the energy density fore- 
most from big bang nucleosynthesis, which is the next 
important milestone in the evolution of the universe af- 
ter the QCD phase transition. Caprini and Durrcr found 
that magnetic fields produced by a causal production 
mechanism (in contrast to magnetic fields produced dur- 
ing primordial inflation) can be strongly limited via their 
inte gra ted energy density and the shape of the spectrum 
[63ll64|. They argued that the spectrum of the gener- 
ated magnetic field must fall off with a steep power law 
for uncorrelated superhorizon scales, i.e. B\ ex A - " with 
n > 2. As stated earlier the typical comoving length scale 
of galactic magnetic fields is 0.1 Mpc which is comparable 
to the shortest magnetic field mode that survives plasma 
damping processies up to recombination {65l l66j . This 
scale is clearly much larger than the comoving horizon 
size H^ 1 ~ 10 pc at the QCD phase transition. There- 
fore even a relatively small field strength at the 0.1 Mpc 
scale requires a magnetic field at the 10 pc scale that is 
larger by many orders of magnitude easily resulting in a 
very large integrated magnetic field energy density. We 
use the bound on an additional radiation e nerg y density 
at big bang nucleosynthesis found by ref. [67[ allowing 
at most 1.6 additional effective neutrino families at the 
98% confidence level. The integrated magnetic energy 
density is thus bounded from above by Bqcd = 5 ■ 10 13 G 
which limits the strength of the comoving seed field to 
.BfTiMpe < 10~ 22 G. Our previous estimate of the gen- 
erated magnetic field consequently does not violate this 
bound, but the field strength is very low and may not 
suffice to seed large scale magnetic fields if not enhanced 
sufficiently. In [681 ] it was found that an inverse cas- 
cade mechanism could transfer some field strength from 
small to larger scales thus partially escaping the effects 



of plasma damping. The inverse cascade mechanism re- 
quires a non-vanishing helicity of the primordial mag- 
netic field, as one can expect in the presented scenario 
due to the large baryon asymmetry, thus one may still to 
successfully seed large scale magnetic fields fields at the 
QCD phase transition. 

VIII. GRAVITATIONAL WAVES 

The final signal of the QCD phase transition that we 
would like to discuss are gravitational waves. The process 
of nucleation and subsequent bubble collisions will stir 
hydrodynamic turbulence producing graviational waves 
in the process [5^, I691 - I72T ] . Again the Hubble parameter 
gives an important scale for the spectrum (5(J |7l| . Since 
the production mechanism is causal a peak frequency has 
to greater or equal to the Hubble frequency, v pea k > v n- 
By how much this peak scale differs will depend on the 
details of the production mechanism and most impor- 
tantly on the relevant time- and lengthscalcs that might 
be significantly different from the Hubble frequency. Let 
us assume an exponential nucleation rate T oc exp(t/r) 
with r being the characteristic timcscalc for the nucle- 
ation process. Then the peak frequency of the spectrum 
due to the collision of bubbles will be given by 

peak \ t J VlSOMeV,/ V 50 J 

(65) 

where T* is the reheating temperature and the result 
is already redshifted to the present day frequency. It 
is common to denote a gravitational wave spectrum 
in terms of a characteristic strain amplitude which is 
definied in the following way 

h c {v) = 0.9 ■ 10- 18 (^) (jfy [LVM] 1/2 (66) 

Using the results of [tT| with the above estimates one 
arrives at a peak strain amplitude for the bubble collision 
peak of 

l/B ^ i-^-is/ T \ 2 /!50MeV\ / 50 V 7 

(67) 

Bubble collisions will also create hydrodynamic turbu- 
lence that will stir gravitational waves with a slightly 
lower peak frequency 

v peak — O- 3 u peak (68) 

but with a higher peak amplitude 

M"£«fc) - 2.1fcc(K£»*) (69) 

for a strongly first order phase transition [7l[ . The addi- 
tion of magneto-hydrodynamic (MHD) turbulence could 
further boost the peak amplitude as shown in [73|, [74j . 
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For frequencies lower than the Hubble frequency the 
spectrum should be uncorrelated white noise. The ap- 
proximate shape of the strain amplitude spectrum is then 
given by 



h c (v) oc v 1 / 2 for v < H 
h c {y) oc v n for v > v. 



peak 



(70) 
(71) 



where the spectral index n should be at least -2 if the 
number of bubble collisions is low. If multi-bubble colli- 
sions play an important role the index is expected to be 
lower where simulations of the collision of vacuum bub- 
bles as done in [72[ find n = —3/2 while older simulations 
found — 2 < n < — 1 [(3jl[7(|. Note that a gravitational 
wave spectrum with n = — 1 would be UV-divergent un- 
less there is a cutoff [82|], so n < —1 for a realistic spec- 
trum. 

The Parks Pulsar Timing Array (PPTA) measures tim- 
ing residuals in pulsar signals to put upper bounds on 
a stochastic gravitational wave background in a rela- 
tively narrow frequency band around 10 -8 Hz (75l [76j . 
The PPTA results already allow to limit the nucleation 
timescale with the presently available data to t/H^ 1 < 
0.12. This limit will improve to r/i?" 1 < 0.06 for the full 
data of the Parkes Pulsar Timing Array project [75| or 
even beyond that depending on the position of the peak 
[77l 178} since the sensitivity depends on the spectral in- 
dex of the gravitational wave background. The results 
can be seen in figures [11] where the expected spectrum 
of gravitational waves for the former case is shown with 
three different high frequency slopes and the approxi- 
mate sensitivity regions of existing and future detectors. 
In [7§] the most optimistic results from [73|, 741 are com- 
pared with the most optimistic results from j69tl70| where 
the former have a slightly higher peak amplitude due 
to MHD turbulence while for the latter the contribution 
from multi-bubble collisions is stronger in the high fre- 
quency regime. 

The planned Square Kilometer Array (SKA) will im- 
prove the sensitivity in ^lg W {u) by about four orders 
of magnitude j8(|. Thus SKA will lower the bound 
on t/H^ 1 by about an order of magnitude as visible 
in figure QTJ If multi-bubble collisions are important 
detection via the spaceborne Laser Interferometer Space 
Antenna (LISA) could also be possible if the high 
frequency spectral index n > —1.4 and r/ff _1 > 10~ 2 . 
In figure [12] the limiting case of a strain amplitude 
spectrum for t/H~ 1 = 0.005 is shown below which the 
signal would be unobscrvablc with cither SKA or LISA. 

Furthermore it has been found that the QCD-phase tran- 
sition will also leave a steplike imprint on the spectrum of 
primordial gravitional waves due to the str ong reduction 
of the radiation degrees of freedom [81| • In [79| this result 
was confirmed also for several lattice equations of state. 
Furthermore the effect of a little inflationary phase on the 
primordial spectrum was examined and a strong power- 
law suppression for frequencies larger than the Hubble 
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FIG. 11: Largest strain amplitude spectrum at t/H^ 1 = 0.12 
that is still compatible with the data of the Parks Pulsar Tim- 
ing Array. A shorter duration of the phase transition reduces 
the amplidude and shifts the peak to higher frequencies. De- 
tection with LISA would only be possible if multi-bubble col- 
lisions play a significant role. 




FIG. 12: Strain amplitude spectrum for t/H 1 = 0.005 below 
which the signal would be unobservable with either SKA or 
LISA. 



frequency at the QCD phase transition was found. 



IX. SUMMARY 

We have reexamined the idea of a little inflation at 
the QCD phase transition and extended the discussion 
as compared to our previous publication fl2l ]. We 
found interesting cosmological implications such as the 
suppression of primordial density fluctuations up to 
dark matter mass scales of M max ~ 10 6 Mq relative 
to the large scale spectrum due to the change of the 
global equation of state. This could have interesting 
consequences for the physics of globular clusters and 
the emergence of the first stars and could also have an 
impact on the cuspy core density distribution of dark 



15 



matter in small galaxies and the too large number of 
halo structures seen in standard structure formation. We 
found that the baryon density can actually be so large 
that one may even expect color superconducting phases 
to be present before the onset of the little inflation, this 
might pose an alternative route of investigation for the 
scenario. 

We also discussed the production of primordial magnetic 
fields that may be strong enough to seed the presently 
observed galactic and extragalactic magnetic fields. 
Furthermore we addressed the production of a spectrum 
of gravitational waves around a peak frequency of 4 • 10~ 8 
Hz that may be observable via pulsar timing in the 
future [ill H(| • Dark matter properties are also strongly 
affected as the annihilation cross section for cold dark 
matter has to be up to nine orders of magnitude lower 
to give the right amount of dark matter today, which 
can be probed at the LHC by detecting the neutralino 
with an unexpected low annihilation cross section. The 
conditions in such a cosmological phase transition would 
then be closer to the situation in heavy ion collisions or 
even the centre of neutron stars than to the standard 
QCD phase transition in the hot big bang scenario. 
Hence, the upcoming FAIR facility would actually for 
the little inflation scenario be a probe for the physics of 
the early universe. 
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X. APPENDIX 

Now let us discuss the analytic solutions to the system 
of differential equations for the perturbations 



Si 



_&^i 1+ *^-3(l + B( )2 (72) 
a Ha a 

- c 2 sl — 8 t - (1 + Wi)— a (73) 

a Ha Ha 

a = - *( 1 + »*> S (74) 
(«) +!(! + ») 

for the most relevant cases. We will limit the discussion 
to stating the approximate equations of motion for the 
radiation and the dark matter component by directly 
giving the dominant analytic solutions. First let us look 
at the growing super-horizon solutions in the case of 



super-horizon: k/H <C 1, 

radiation domination: <? s — w — 1/3, 

a — —Sr/2, k/H oc a 


S R - ^R + 1 & R 


S'dm - h 5 R 

*{>'dm - + km 5 * 


&* = 4(&) S 


Sdm = §<5r 
4>dm = %i>n 



TABLE II: Equations of motion and dominant perturbation 
solutions for a radiation dominated universe in the super- 
horizon limit. The constant A fixes the amplitude at hori- 
zon crossing; and is scale independent for a Harrison-Zeldovic 
spectrum [47L I4H ] . 



a radiation dominated universe, which sets the initial 
condition for the numerical calculations. These solutions 
set the relevant initial conditions for our numerical 
calculations. The approximate equations of motion and 
the resulting solutions are summarized in table iHl 

Inflationary solutions 

Now we want to additionally find the solutions for the in- 
flationary phase in the different spectral regimes. For the 
inflationary regime it will be most important to examine 
the case of q = 3 because radiation will be less abundant 
than matter soon after the onset of inflation for relevant 
inflation lengths. First of all we need the mean quantities 



1 



w 



'DM 



& 



1-31 w~-l 
. (75) 

and we also need to remember that k/% oc 1/a in the 
following. Now let us examine the two most relevant 
spectral regimes namely the intermediate and the unaf- 
fected regime. For the spectral range we have examined 
numerically none of the modes will stay sub-Hubble 
sufficiently long during inflation to approach an analytic 
limit. This would only be the case for modes that stay 
similar or even below the Hubble frequency for the whole 
duration of inflation. It turns out that the solutions 
in this case are combinations of Bessel functions that 
cannot be found by simple analytic means. Thus we 
skip a lengthy discussion for these modes and directly 
jump to the other two regimes that are relevant to 
the discussion and have analytic solutions that can be 
derived rather quickly. 

First let us discuss the intermediate modes with 
1 > I > (1 I w) 1 / 2 . The approximate equations of 
motion in this regime and the corresponding solutions 
are summarized in table IIIII The solutions for the 
dark matter perturbations imply that 6 dm and 5 a will 
be frozen very quickly until the end of inflation and 
approach a constant value. Note that (1 + w) 1 / 2 drops 
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intermediate modes: 


1 > i > (l + w) 1/2 , 


inflationary phase: c 


2 ~ 1/3, it; ~ —1, k/H oc 1/a, 


«--3(£) 2 (^) 3 * 


DM 


&R — 3$DM 


S'dm ^ 9 (f ) 2 %5dm 


,/)' ~ 1 k A' 
™R — 3H°DM 




8r = |5dm 


fcM=C iex p [- 9 (f) 2 (^) 3 ] 







TABLE III: Equations of motion and dominant perturbation 
solutions for intermediate modes during the short inflationary 
phase, Ci and C2 are constants. 



quicker than ^ so any mode that enters this regime 
stays there until the end of inflation. 

Now let us turn to the unaffected modes that are given 
by the condition 1 3> (1 + w) 1 / 2 3> jr. The correspond- 
ing approximate equations of motion and the resulting 
solutions in this regime are again summarized in table 
IIVI Note that the solutions for 5r and 8dm are exactly 



unaffected modes: 1 » (1 + w) 1/2 » J|, 


inflationary phase: c 2 ci 


i 1/3, w ~ — 1, fc/W oc 1/a, 


OL ~ — |(5DAf 




R ~ 3° DM 


8 DM — ffojlf 


?,/ Iks ,8k 

VR - - affB + 9to 


feM — -\"ij>DM + f^foM 


Sr = ^8 DM 


odm = Ba 2 cx (f ) 2 


IpR = |?/>DM 


^M = |iBa 2 oc(f) 



TABLE IV: Equations of motion and dominant perturbation 
solutions for unaffected modes during the short inflationary 
phase, B is a constant. 



the same solution as for the radiation dominated super- 
horizon case as found before. Comparing these results 
to the analytic super-horizon solutions in the radiation 
dominated universe we find that 8 dm and Sr have the 
same growing mode oc a 2 , thus the naming of the spec- 
tral region as "unaffected" is justified. On the other hand 
ipDM and if; n grow only linearly with the scale parameter 
in contrast to a cubic growth in the radiation dominated 
super-Hubble case. This is actually necessary to keep the 
spectrum scale invariant on large scales (48j since this re- 
quires Si/ipi oc k/T-L. The latter keeps the amplitude at 
horizon entry independent of the wavenumber for large 
scales. 

Furthermore the comparison between the intermediate 
and the unaffected cases shows that one may expect a rel- 
ative suppression in the fluctuations below the unaffected 
part of the spectrum because the intermediate modes are 
frozen in while the unaffected modes grow in the respec- 
tive limit. 
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